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Abstract 
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1 Introduction 



Recently the inclusive decay spectra of heavy hadrons have attracted renewed 
attention. It has been shown that by a combination of Heavy Quark Effective 
Theory and the Operator Product Expansion one may perform a systematic 
l/mg expansion of the charged lepton spectrum in inclusive heavy hadron 
decays - [@] • However, it turns out that an adequate description of the end- 
point region, where the charged lepton energy Ee ~ Emax, requires a partial 
resummation of the Operator Product Expansion, yielding a result analogous 
to the leading twist terms in deep inelastic scattering. In particular, it in- 
volves the analogue of the parton distribution function, which in the present 
case describes the endpoint region of the lepton spectrum. While in most of 
the phase space the nonperturbative corrections may be described in terms 
of a few parameters, the endpoint region needs the input of a nonperturba- 



tive function @|-[|ri|]. These recent ideas constitute a substantial conceptual 
progress towards a model independent description of inclusive decay spectra, 
including the endpoint region. 

However, aside from the nonperturbative effects one also has to include 
QCD radiative corrections before one may confront the theory with data. 



These corrections have been calculated IIT2|1-|[TO1 and are known to be small for 



b ^ c transitions over the whole phase space. This is also true for the b ^ u 
case, as long as one is not too close to the endpoint E^ ~ Emax- Neglecting 
the mass of the u quark one finds Sudakov-like double logarithms of the 
form In {Emax — Ei) as well as single logarithms hi{Emax — Ei), indicating a 
breakdown of perturbation theory close to the endpoint. 

In the present context the QCD radiative corrections have been studied 
by Falk et al.|T^ and Bigi et al.|TO[. Falk et al. arrive at the conclusion 
that the endpoint region of the 6 — > m decays is strongly affected by QCD 
radiative corrections and in principle a resummation to all orders becomes 
necessary. The conclusion of [1^ is that the extraction of the nonperturbative 
effects becomes practically impossible due to large and mainly uncontrollable 
perturbative effects. 

In the Sudakov logarithms are exponentiated leading to a strong 
modification of the endpoint region. The authors of suggest that un- 
controllable radiative corrections discussed in may be eliminated by the 
appropriate choice of the scale entering the strong coupling. They use 
/i ~ Arrih, where A is the characteristic momentum of light degrees of free- 
dom in the heavy meson. However, this leads to a strong modification of the 
shape of the spectrum near the endpoint and the reliability of the calculation 
becomes questionable. 

Unfortunately, it is only the endpoint region of charmless B decays, which 
is not buried in the huge background of the charmed decays. The kinematic 
endpoint of the charmed B decays is (m^ — m|,) / {2mB) ~ 2.3 GeV, while for 
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the charmless semileptonic decay it is (m^ — m^)/(2mB) ~ 2.6 GeV. Hence 
there is only a window of at maximum 300 MeV, which may be attributed 
solely to the transition b — >■ uiu. The size of this window is of the order of 
the mass difference of the heavy meson and the heavy quark, A = ms — rub, 
and thus it is not dominated by a few resonances. This would be the case in 
a smaller region of the size A(A/mf,), which is of the order of tens of MeV, 
and in which methods mentioned above would fail, since they are based on 
parton hadron duality. 

The present paper is an attempt to define quantities which are on one 
hand experimentally accessible, and on the other hand allow to disentangle 
perturbative and nonperturbative effects. Our suggestion is to calculate ap- 
propriate moments 7W„ of the measured spectrum, which are defined in such 
a way that they allow an expansion in both ^^(mf,) and A/rrib, at least for 
some range of indices n. 



Similar ideas have been put forward previously in and [|10]- However, 



in these papers there has been no detailed numerical study of such moments, 
including both radiative and nonperturbative effects. 

In the next section we shall reconsider the perturbative and nonperturba- 
tive contributions to the decay spectrum of inclusive semileptonic charmless 
B decays and give the definition of the moments. In section 3 we perform 
a numerical study from which we obtain constraints on the range of the in- 
dex n of the moments. Finally we give our conclusions and comment on the 
extraction of Vub- 



2 Perturbative and Nonperturbative Contri- 
butions to the Spectrum 

It has been shown in -[§ that one may obtain a systematic 1/mb expansion 
for inclusive decay rates and the corresponding spectra. The leading term in 
this expansion is the free quark decay, and the first nonvanishing corrections 
appear formally at the order These are genuinely nonperturbative 

corrections, which are given in terms of two matrix elements 

{B{v)\K{iDfK\B{v)) = 2mB\i (1) 
{B{v)\K{-i)a^SDn{^D^)hv\B{v)) = emeAa. (2) 

The parameter A2 is given in terms of the 0~ — 1~ mass splitting A2 = 
(mg*— m^)/4 ~ 0.12 GeV^. The matrix element Ai is not as easily accessible, 
in particular it has not yet been determined experimentally. The theoretical 
estimates vary over a broad range ||18|, |19| of about Ai ~ —(0.3 — 0.6) GeV^; 
thus we shall consider below a typical "small" value Ai = —0.3 GeV^ and a 
"large" one, Ai = -0.6 GeV^. 
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The result for the spectrum of the inclusive decay B 
the mass of the u quark, is given by 
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mb 1927r3 
The 6 function singularities at the parton model endpoint y = 1 indicate a 
breakdown of the operator product expansion. In fact, it has been pointed 
out in H, H that for the decay spectra the expansion parameter is in fact 
A/[mb(l — y)], which becomes large close to the endpoint. Still these terms 
have an interpretation as being the first few terms of a moment expansion of 
a nonperturbative function describing the behavior of the spectrum close to 
the endpoint M-n. 
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where 5^^^^ denotes the fcth derivative of the 5 function. Note that ai vanishes 
and 02 may be read off from (H) to be 02 = — Ai/3. 
The function F{x) may be written formally as 



Fix) 



2m I 



■{B{v)\hJ[x- 



K\B{v)), 



(6) 



where is the positive light cone component of the covariant derivative. 
This function has support only in a region of x ~ 0(1). In particular, the 
function F leads to some "smearing" of the endpoint region, with the effect 
that the endpoint of the spectrum is shifted from the parton model endpoint 
mb/2 to the physical endpoint ms/2, where is the S-meson mass. Below 
we shall use some model input for F to estimate the size of the higher order 
corrections close to the endpoint. 

Aside from these nonperturbative corrections there are also perturbative 
ones, which have been calculated some time ago. The order as corrections 
are known for all values of the lepton energy |]T2[-|T^; however, since we 
are interested only in the behavior close to the endpoint, we shall consider 
here only the contributions relevant in this region. It turns out, that the 
order as corrections exhibit doubly and singly logarithmic divergencies at 
the endpoint. Up to terms vanishing at the endpoint, the result reads |]16| 

\v?{l-y) + — \Yi{l-y) + - + n- 
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where at the endpoint 



1 cir(o) 
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The scale fi of the strong couphng is not yet fixed in this expression; 
any scale choice will formally only affect subleading terms. However, from 



physical considerations one may be lead to chose /i 



mt 



y) as it was 



done in [T^Q, leading to an even more singular behavior as ?/ ^ 1. In any 
case, a resummation of the singular terms becomes mandatory in order to 
describe the spectrum close to the endpoint. 

The common folklore is that the doubly logarithmic terms in (|^ expo- 
nentiate, and at the level of the Sudakov logarithms the radiative corrections 
become 
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Again the problem arises of the scale of and we shall compare below two 
choices. Using the one loop expression for as 
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H = rrih and /i ~ 'mlil — y), the 
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In the endpoint region one obtains qualitatively a behavior of the form |T0 

dv c/r(o) 

dy dy 

where 

. 8 
Co - 1 = 



(12) 



ln(l - y) 



25\n[iml/AlcD)a-y)] 

behaves approximately like a constant in the region of interest |]TD||. How- 
ever, after such a resummation to all orders of perturbation theory, it is not 
obvious, how to disentangle perturbative and nonperturbative contributions 
any more. This becomes clear, if one really treats eo as a constant, in which 
case one may rewrite (|12D 
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^ Actually, the scale chosen in [|iO| is fi^ ^ 
since in the endpoint region I — y ^ A/mb. 



Arrib, but this is equivalent to ^ ^ m^(l — y), 
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which in this form looks hke a nonperturbative contribution. 

The above discussion shows that the endpoint region is in fact difficult 
to describe. On one hand there are large nonperturbative contributions, 
on the other hand there are perturbative contributions, becoming large in 
the endpoint region. These remarks, however, apply only to the spectrum 
itself. It has already been shown in ||^, |10| that it may be helpful to consider 



appropriately defined moments of the spectrum; e.g. the moments of the 
spectrum taken with respect to the parton model endpoint are related to the 
coefficients of the most singular 6 function like singularities at the endpoint, 
in each order in the l/rrib expansion. 

Here we propose to consider a slightly different set of moments, namely 

l+A/m6 

dy 2/"^, (14) 

^ 

which may be rewritten in terms of the moments as defined in 0. As we 
shall discuss in detail below, these moments have, for some range of n, a 
simultaneous expansion in as{mb) and A/rrib. 

This range of n is specified by two requirements. Experimental informa- 
tion will be available only close in the endpoint region, and thus n has to be 
large enough to be sensitive to this region. On the other hand, the larger n is, 
the stronger the sensitivity to the details of the endpoint region is. In other 
words, only for not too large n one may perform the simultaneous perturba- 
tive and nonperturbative expansion; the problems present in the spectrum 
will reappear in the behavior of the moments for large n. 

In fact, from order of magnitude considerations one would expect that 
there is no such range in n, where n is large enough to be sensitive to the 
endpoint, and in which nonperturbative and radiative effects still remain 
under control. However, it turns out from the more detailed study presented 
below that there might be such a a window in n. 

In the next section we shall perform a numerical study for the radiative 
corrections and the nonperturbative contributions, in order to find a range 
in n where one may reliably calculate the moments A^„. 

3 Numerical Discussion 

We shall first consider the moments of the differential distribution in the 
naive parton model, without any corrections. The spectrum is given by the 
parton model expression and one obtains for the moments to the zeroth order 

\A0) 2(n + 6) 
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However, we are mainly interested in the endpoint region, where the parton 
model spectrum may be approximated by (§), and the moments obtained in 
this approximation are 

A^f^r.^ (16) 

In Fig.|l| we plot the moments obtained from the full parton model rate versus 
the approximation (plGl). The comparison between the two sets gives some 
impression, at what n one is sensitive only to the endpoint region; from the 
figure one reads off that already at n ~ 4 one mainly obtains information on 
the endpoint. 

On the other hand, data will be available in the near future only for lepton 
energies above 2.3 GeV. Thus we shall in the following also consider moments, 
in which the integration over n is restricted to a range Uq < y <1 + K/ 

/ai 
dy y"^, (17) 

where realistic values for y^ will be in the region of yo = 0.9. Introducing 
such a lower cut will change (^) to 

= (i - yr) , (18) 

from which we estimate that for y^ = 0.9 the sixth moment will get about 
52%, the tenth moment already about 70% contribution from the region 
yo < y < 1, at least in the naive parton model. On the other hand, an 
upper limit of n is given by the experimental resolution. If data on 6 — m 
semileptonic transitions is available only in the small window between 2.3 
GeV and 2.6 GeV one may not expect enough data to extract moments higher 
than about ten. In what follows we shall thus concentrate on moments less 
than ten. 

Next we turn to the nonperturbative corrections. They have been been 
given in (|^) to order I/ttlq. Taking the moments of one obtains in the 
endpoint approximation (^ 

Mn = Ai<»' + r. 



3m? Vn + 1 J ' ml \n + 1 



(19) 



The dependence on n of the various terms in Ain reflects, how singular the 
contribution is at the endpoint of the lepton spectrum. The most singular 
term is the one with the derivative of the S function; this leads to a linear 
dependence in n. The terms behaving like a S function yield constant terms 
in the moments, while terms with step functions will decrease as 1/n if n 
becomes large. 
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Figure 1: Comparison of the parton model moments for the 
full y dependence (full dots) versus the endpoint approxima- 
tion (open circles). 



In a similar way as for the parton model we may estimate the effects of 
a lower cut off according to ([TtD . For the terms involving step function 
the result is qualitatively the same as for the parton model, while for the 
contributions with 5 function like singularities there is no dependence on the 
lower cut, since these are concentrated dX y = 1. 

Higher moments will become more and more sensitive to what happens 
in the endpoint region. Thus (|TU|) is not valid for n becoming too large. In 
order to estimate the value of n, where the expansion breaks down, we have 
to consider an even more singular contribution. This is contained in the next 
order of the l/mq expansion, and has the form 

where the parameter as is related to the matrix element 
{B{v)\K{tD^){tvD){tD'')K\B{v)). 
Taking the moment of this gives a contribution quadratic in n 

Mnr^Tbn{n-l)^. (21) 



7 



Neglecting the higher orders in the l/mq expansion of the moments is only 
justified, if 



n[n — 1) — ^ < n — 5-. 

mt mi 



Estimating 03 to be (lAil)^/^ _ (590 MeV)^ one obtains 



n < ~ 10, (22) 



which means that only the nonperturbative contributions to the first few 
moments may be calculated to lowest order in the l/mq expansion. It is a 
fortunate circumstance that this upper limit in n coincides with the maximal 
n accessible in experiment. 

In order to estimate higher moments one has to perform a resummation of 
the most singular terms at the endpoint to all orders in I/itlq, in other words, 
one has to include these effects using the nonperturbative function F{x) 
defined in (|^). In terms of this the rate is given by (^, which describes the 
behavior of the lepton spectrum close to the endpoint, namely over a region 



of the order 1 — y ~ y |Ai|/mf,. There is, however, an even smaller region 

I — y (y^lAil/mb)^, the resonance region, in which only a few resonances 
contribute to the spectrum. This region will start to contribute significantly 
to the moments for 

n ~ (^L)2 _ 100 (23) 

which means that moments with n > 100 will be determines from the con- 
tributions of only very few light resonances. 

The function F is genuinely nonperturbative and it has been considered 
using several models. One frequently used model is the one of Altarelli et 
al. |T^, which has been rewritten in terms of the nonperturbative function 



F(x) ||2^ given by 




F(x) = ^^exp|-^(-^l^-^(l-x)) ). (24) 

Here A = ms — rrib, where is the mass of the 5-meson. The second 
parameter is the so called Fermi-momentum which corresponds to the 
motion of the heavy quark inside of the heavy meson; see |]2D| for precise 
definition. Finally, p is implicitly given in terms of the other two parameters 
by 

-^^P^P^"JM2)^ (25) 
where Ki is a modified Bessel function. 
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The model as given in has only one parameter p, which has a physical 
interpretation, namely p = ""^ipectator/PF- may also relate the model 

parameters to the QCD matrix element Ai, since appropriate moments of 
F are related to matrix elements of powers of covariant derivatives between 
heavy meson states. In this way one obtains the relation 

^Tiikr' '''' 

which may be used to determine the value of the model parameter p in 
terms of A and Ai; it has been pointed out in pO| that only the combination 



^ = — Ai/(3A^) enters in the model. In particular, eq.(PBD has no solution, 
if ^ is larger than 0.57. For a value of A of 500 MeV, this means that the 
model cannot accommodate a value for — Ai of more than 0.42 GeV^. For 
this value we have p = 0, and the nonperturbative function simplifies 

F(x) = -exp ("- (^-^^' V (27) 
71 y n J 

We shall use this model to estimate the effects of the most singular terms 
occurring in higher orders in the l/rrib expansion. 

The value of Ai is not yet known accurately. We shall consider below 
a "small" value Ai = 0.3 GeV^ as well as the largest value possible in the 
above model, Ai = 0.43 GeV^. The first value of Ai corresponds to p = 0.35, 
while the large one is the limiting case p = 0. In our final result we shall also 
consider a "large" value Ai = — 0.6GeV^. Furthermore, we will use rUb = 4.7 
GeV and A ~ 500 MeV in the following numerical analysis. 

In Fig.|^we plot the nonperturbative corrections 5Ai^^^^ to the moments, 
§M'^P) = _A/j{np) _ _A^(o). The upper plot corresponds to Ai = -0.3 GeV^ 
and the lower one to the maximal value for Ai which can be accommodated 
in the model (24), Ai ~ —0.43 GeV^. The solid dots are the nonperturbative 



corrections according to (|T9|), while the solid triangles are the corrections 
using the model (|2^) for the nonperturbative function F. 

We have also plotted separately the contributions to the moments origi- 
nating from the term linear in n of (|19D (open circles) and the rest, i.e. the 
constant terms and terms decreasing with n (open squares). 

In both cases, small and large Ai, one finds a substantial contribution 
from the chromomagnetic moment term A2 at least for small moments. The 
term from the kinetic energy operator has a linear dependence on n (c.f. 
(|19D) while the chromomagnetic term behaves like a constant for large n. The 
linear terms of ( ^^) are well reproduced by the model for the nonperturbative 
function; this is to be expected, since the nonperturbative function contains 
only the most singular contribution in each order of the l/rrib expansion, i.e. 
the leading twist term. On the other hand, from the fact that the linear term 
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Figure 2: The nonperturbative corrections to the moments 
for Ai = -0.3 GeV2 (upper figure) and Ai = -0.43 GeV^ 
(lower figure). The sohd dots are the corrections using the 
lowest nontrivial contributions (|1^), the solid triangles are 
the corrections from the nonperturbative function as given in 
(|27|), the open circles are the contributions from the nAi alone 
in ([19|) , and the open boxes are the contributions except the 
nAi term. 
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Figure 3: The nonperturbative corrections to the moments, 
including a lower cut of uq in the integration as in ([T7|). We 
use p = 0, corresponding to Ai = —0.43 GeV^. 

in n from (|19D already approximates the nonperturbative function quite well, 
one may conclude that the expansion ([T9|) is in fact sufficient for the range of 
n we are considering. In other words, the leading twist contribution modelled 
by the ansatz (p^) may be replaced by the term proportional to Xi/ml in 
the l/ruf, expansion of the moments with high accuracy. 

However, the result (^) differs from what the nonperturbative function 
gives, mainly because of the chromomagnetic moment term and the constant 
terms proportional to Ai. These term are less singular than the contribution 
of nAi, but they contribute substantially to the moments n < 10. On the 
other hand, the leading twist terms indicate that the l/m-b expansion is 
satisfactory for n < 10 and we conclude the expansion (|19D is justified for 
the moments in the range considered. 

Finally, one may also consider the effect of a lower cut off i/q in the 
integration over y as in (|T7|) . The strongest effects are expected for large Ai, 
so we plot in fig.^ the nonperturbative corrections for p = in the model 
(plD , for different values of the cut off yo. For the realistic case yo = 0.9 one 
finds only small corrections for 4 < n < 10, and the conclusion concerning 
the minimal n obtained from the parton model remain valid. 

Next we consider the radiative corrections. In Fig.^ we plot these correc- 



11 



_Q ^ I \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
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n 

Figure 4: The radiative corrections to the moments. Open 
circles correspond to (|^) with as{mb) = 0.26, the crosses cor- 
respond to (H) with as{mb) = 0.26, and the x'es are obtained 
from ([n]). 

tions to the moments 6A4^^'^^ obtained from the one loop result (|^), from the 
exponentiation of the leading double logarithms, for both fixed scale /i = rub 
and for variable scale = mb{l — x). We use in (0) the value of as at rub, 
ctsimb) = 0.26, corresponding to Aqcd = 250 MeV. 

The radiative correction contribution to the moments with > 4 is rather 
stable, in other words, it does not depend on how one deals with the endpoint 
singularities. Compared to the nonperturbative corrections for small Ai they 
are much larger, dominating the corrections to the moments. This smallness 
of the nonperturbative corrections is due to an almost cancellation between 
the contributions of Ai and A2 in the region of n we are considering. 

The effect of a lower cut in the y integration for the moments as in (p!7|) 
leads to a doubly logarithmic dependence of the form ln^(l — ?/o) on the lower 
cut Ho, and, if yo comes close to 1, the effects of the lower cut may become 
large, forcing us to go to higher n. However, numerically it turns out that 
with realistic values of yo the situation still remains under control. This is 
shown quantitatively in fig.^, where we use expressions (upper figure) 
and d^) (lower figure) to estimate the effect of a lower bound. Again, for the 
realistic value yo = 0.9 the corrections to the moments 4 < n < 10 are less 
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Table 1: Numerical values for the functions A(n) and a(n) 



than ten per cent. 

Our main conclusion is thus that one may use the combined expansion 
for the moments below n = 10, which is given by 
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The values for A(n) and (j{n) are tabulated in tab.|l]. 

Furthermore, the fact that data will be available only for y > yo 0.9 only 
has a small effect on the moments in the range of n considered; the moments 
including a lower cut yo may still be treated in the combined expansion. 

The final result for the moments is plotted in Fig.^, where the moments 
obtained from the combined expansion in as{mb) and A/nib according to (^) 
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Figure 5: The effect of a lower cut off yo (c.f. (|T^) on the 
radiative corrections to the moments. The upper figure are 
the moments obtained from (|^) with lower cut off, while the 
lower plot is obtained from (||). 
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Figure 6: The final result for the moments, including both 
radiative and nonperturbative corrections according to (pHf). 
The solid dots are the result for Ai = —0.3 GeV^, the solid 
boxes are for Ai = —0.6 GeV^. For comparison we also plot 
the parton model result (open circles). 

are shown. The solid dots are the result for "small" Ai = —0.3 GeV^, and 
the solid boxes are for "large" Ai = —0.6 GeV^. 

To summarize, one may analyze the moments up to 12 before the cor- 
rections become too big, of order 100%. We also see that nonperturbative 
corrections are very small for small values of Ai, and key role is played by 
the radiative ones. The reason for the smallness of nonperturbative correc- 
tions is an almost complete cancellation between the terms proportional to 
Ai and A2. For the case Ai = —0.6 GeV^ the cancellation is only partial 
and nonperturbative correction are sizable and positive. This indicates that 
higher twist effects may play an important role also in the endpoint region, 
and taking into account only the leading twist contribution, corresponding 
to the function F, is not enough. This remark, however, applies to a precise 
description of the spectrum in the endpoint region, while for the moments 
with sufficiently small n it is safe to use the combined expansion in as{mi,) 
and A/rrih, which contains pieces of nonleading twist in the constant terms 
and the contributions decreasing with n. 
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4 Conclusions 



Since data on 6 — m transitions will be restricted to a small window between 
2.3 and 2.6 GeV, it is mainly the endpoint of the spectrum which will be 
accessible in experiment. On the other hand, the endpoint region is the most 
difficult region from the theoretical point of view, since here large radiative 
corrections are entangled with large nonperturbative effects. In fact, it is 
not even obvious, whether and how one may distinguish the two sources of 
corrections close to the endpoint. 

The main result of this paper is that one may define suitable averages 
of the inclusive decay distribution of semileptonic b u transitions, which 
on one hand may be calculated reliably, and which on the other hand are 
mainly sensitive in the window between the kinematic endpoints of 6 — > c 
and b ^ u semileptonic decays. We have concentrated on moments of the 
energy spectrum of the charged lepton and have performed a detailed numer- 
ical analysis, which shows that moments for ri < 10 may be systematically 
calculated in a combined as{mb) and A/rrib expansion; both types of correc- 
tions may be studied systematically. Furthermore, the moments with n > 4 
are sensitive mainly in the experimentally accessible window. Thus there is 
indeed a region 4 < n < 10 for which the moments defined above may be 
useful to perform a determination of \Vub\- 

In fact, such a combined expansion is even valid for moments defined 
with a lower cut off yo, as given in ([T7| ) as long as yo is not too close to one. 
Furthermore, for a realistic value yo = 0.9 the deviations of the moments 
including a cut from the full ones are small in the region of n considered. 
Thus for an experimental analysis one may as well deal with the moments 
including a cut. 

In order to extract \Vub\ from these moments one has to know the two 
matrix elements Ai and A2. While A2 is known from the hadron spectrum, Ai 
is still uncertain; we have used typical "large" and "small" values, but one 
may hope to extract this parameter from experiment eventually. 

One possibility to perform a determination of Ai is in fact to use the 
moments defined above. The ratio of two moments will depend only on 
Ai/mf, X2/rnl and as{mb). The strong coupling and A2 are known, and 
one may in principle extract Ai and rrib from two ratios of moments. After 
having done this one may then use any one of the moments to obtain and 
hence \Vub\- To proceed along these lines one needs a measurement of three 
moments in the range 4 < n < 10; any measurement of additional moments 
may be used to cross check the n dependence of the moments. 

The idea presented here may be refined in various aspects, and we consider 
this paper as a ffist try to study whether one may extract \Vub\ using the 
moments of the decay distributions. As far as the radiative corrections are 
concerned, we have only considered the singular and nonvanishing terms 
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close to the endpoint. In this point one may refine the analysis by taking 
into account the complete corrections, which may be found in the literature. 



Furthermore, also the corrections of order A/rrif, have been considered [21 
and may be also included in the analysis. However, we do not expect that 
any of these higher order corrections will change any of our conclusions. 

Finally, one may also think of considering other averages of decay distri- 
butions and generalize the moments according to 

Mn= J dyj^niy)j- (31) 

^ 

where JF„ is some set of functions. A good choice of JF„ may possibly allow 
a more reliable calculation of the corresponding moments Ain , and on he 
other hand be more sensitive to the experimentally accessible region. 
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